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Abstract 

Autocorrelation of two pions produced in heavy-ion collisions at intermediate pt 
is calculated in the framework of the recombination model. The differences of the 
pseudo-rapidities and azimuthal angles of the two pions are related to the angle between 
two shower partons in a jet. It is shown how the autocorrelation distribution reveals 
the properties of jet cone of the shower partons created by high-p^ partons in hard 
collisions. 

PACS numbers: 25.75. Gz 



1 Introduction 

Correlations among hadrons produced at high and intermediate p? in heavy-ion collisions 
(HIC) have generated considerable interest in their implications on how the hadronization 
mechanism and jet structure may differ from those in pp collisions [l]-[4]. Experimental 
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investigations in the subjects can broadly be divided into two types: those that use triggers to 
identify near- and away-side jets [1, 2, 4, 5] and those that use no triggers [3, 4, 6, 7]. Among 
the latter is the study of autocorrelations. Autocorrelation has been used extensively in time- 
series analysis. Its application to multiparticle production in HIC was pioneered by Trainor 
and his collaborators [3, 4, 8, 9], and has generated a wealth of information independent of 
triggers and their biases [6, 7]. Theoretical interpretation of autocorrelation has been slow in 
its development. In this paper we present the first prediction of how autocorrelation should 
behave at intermediate px in the framework of parton recombination [10]. 

Theoretical studies of hadron correlation follow a wide variety of approaches; they differ 
mainly in the ways in which the interaction of jets with the medium is treated. Energy 
loss of hard partons in medium has been extensively studied in pQCD [11, 12, 13, 14, 15], 
and its application to correlation on the same side of a trigger particle has been considered 
with emphases placed on various separate but related issues, such as the medium- modified 
fragmentation functions, the effect of collective flow on jets, soft hadrons associated with 
medium-induced radiation, etc. [16, 17, 18, 19, 20, 21]. In such studies hadronization is 
usually treated by use of local parton-hadron duality (LPHD) or in terms of fragmentation 
functions (FFs), neither of which can shed any light on the baryon puzzle in the intermediate 
Pt region, e.g., the high p/n ratio [22]. Since the baryon problem has been well explained 
in the parton recombination model [10, 23, 24, 25], the approach that we adopt here for the 
study of hadron correlations will be in the framework of that model. In that framework a 
number of investigations have already been carried out on dihadron correlations [26, 27, 28, 
29]. The medium effect is taken into account through the recombination of shower partons 
with the thermal partons that are in the environment, when the hard collision occurs near 
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the surface of the medium. Our emphasis in this paper is on autocorrelation. Correlation of 
the hadrons produced on the side opposite to the trigger particle has also been investigated 
recently in numerous ways [30, 31, 32, 33, 34]. Since the path length of a hard parton in the 
medium in such problems is large, they belong to a category of subjects not immediately 
relevant to our concern here, which is restricted to hadrons that are in the same jet. 

The use of trigger has its advantages, especially in showing the properties of the away- 
side azimuthal distribution that reveals the effects of jet quenching. However, it is necessary 
to subtract the background, which is not unambiguous. Autocorrelation is a measure of the 
difference between two nearby values of a variable, with all other variables being integrated 
over. When those values are close, they are dominated by contributions arising from the 
same jet in an event. No background subtraction is needed. So far autocorrelation in the 
data from the Relativistic Heavy-Ion Collider (RHIC) has been analyzed for differences in 
pseudorapidity rj and in azimuthal angle 0, but only at low px [6, 7]. The model that we shall 
use to study the autocorrelation in jets involve thermal-shower recombination for which the 
reliable px region is above 2 GeV/c. Thus at this point our predictions cannot be compared 
to the results of the autocorrelation analysis of the experimental data. Nevertheless, on 
theoretical grounds it is of interest to show how the angular distribution of shower partons 
can be related to the autocorrelation of pions in the differences in 77 and 0. We await with 
anticipation the relevant data that are forthcoming. 
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2 The Problem 

Autocorrelation is a measure of two-particle correlation in multiparticle production with a 
minimum loss of information and without trigger bias. It has been studied extensively in pp 
and HIC [6, 7, 8]. At intermediate to high px the dominant contribution to the two-particle 
correlation comes from jets, when the momenta for the two particles are close together. Our 
problem in this paper is to calculate the autocorrelation distribution from a model in which 
the shower partons in a jet has certain prescribed properties, and in this section we outline 
our plan of attack, leaving the details to the following sections. 

Let X; L be an attribute of the momentum pi of the i th particle, such as its angle relative 
to some axis. The two-particle correlation function in terms of Xi is 

C 2 (x 1 ,x 2 ) = p 2 {x u x 2 ) - pi(x 1 )pi(x 2 ) (1) 

where p\ and p 2 are one- and two-particle distributions, respectively. For autocorrelation 
one defines the sum and difference of Xi 

x± = x 2 ± xi (2) 

and rewrite Eq. (1) in the form 

C 2 (x + ,x_) = p 2 (x + ,x_) -pi((x + -x-)/2)p 1 ((x+ + x-)/2) . (3) 

This would vanish if there is no correlation in p 2 . Anticipating correlation in the variable 
x_ and mild dependence on x + , one defines the autocorrelation distribution to be 

Mx-) = 4 / dx+C 2 (x+,X-) , (4) 
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where the integration is carried out over a range R. Usually, if there are other variables in the 
problem, which are, however, not germane to the correlation measure, they are integrated 
over also. If the range R is wide enough so that the boundaries depend on X-, then one 
must proceed carefully to account for the x_ dependence arising from R. Details of that 
problem related to binning can be found in [8, 35]. In our consideration in the following 
R will be small enough not to involve such complications. Clearly, only the correlated part 
in p 2 (x + ,xJ) contributes to A(xJ). If p 2 (x + ,xJ) has mild dependence on x + , not much 
information is lost by the integration in Eq. (4). A(xJ) treats the two particles on equal 
footing and requires no subtraction of background besides what is explicit in Eq. (3). 

In HIC if the transverse components of pi are > 2 GeV/c, then jets are involved; further- 
more, if the angle between pf and p 2 is less than, say, 7r/4, then the two particles are highly 
likely to be the particles in the same jet. Thus the angular differences between the momen- 
tum vectors provide information about the structure of the jet. Let the angular variables of 
Pi and p2, referred to the longitudinal axis, be (#i,0i) and (9 2 ,<f> 2 ), respectively. Define 



In a central collision + is irrelevant by azimuthal symmetry, which we shall assume. Thus the 
essential variables for correlation are 9 + , Q- and 0_. If the correlation function is determined 
experimentally at midrapidity with a narrow rapidity window, the range of 9 + is not large. 
Then in applying Eq. (4) to this problem, we have 



where Re + is some range in 6 + that is not constrained by the ranges of 6L and 0_ of interest. 
Experimentally, the pt variables are integrated over specific ranges of choice, and are not 



6± = 6 2 ± 6 



<t>± = 02 ± 01- 



(5) 




(6) 
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expressed explicitly. It is important to note that A(6L,0„) depends only on the difference 
angular variables, and is therefore, in principle, independent of the coordinate system in 
which the angles are defined. The only angle associated with two momentum vectors that is 
independent of the coordinate system is the angle x between the two vectors, i.e., 

cos x = cos Q\ cos 6> 2 + sin 9\ sin 6 2 cos (0 2 — 0i ) , (7) 

which can be expressed in terms of 9± and 0_ as 

cosx = \ [cos#_(l + COS0_) + cos# + (l — COS0_)] . (8) 

How this is used in Eq. (6) to obtain the autocorrelation distribution will be discussed in 
Sec. 4. 

The angle x provides the crucial link between the observables and the variables that can 
suitably describe the dynamics of the partons that are associated with jets. At intermediate 
and high p T in HIC the partons that hadronize are the shower partons in jets produced by 
hard scattering. The axis in reference to which those partons can best be described is the jet 
axis. For detection at midrapidity the jet directions are approximately perpendicular to the 
beam axis. Thus to relate the angular variables of the partons, referred to the jet axis, to 
the angular variables of the hadrons, referred to the beam axis, is a complicated geometrical 
problem. In the recombination model the hadrons and the constituent partons are collinear. 
The angle between two shower partons is therefore also the angle between the two pions 
that they hadronize into. That angle is x- Hence, x serves as the bridge that connects the 
momentum space of the partons in the underlying dynamics to the momentum space of the 
hadrons that can be measured. Autocorrelation in terms of 9^ and can thus directly 
reveal the angular properties of the shower partons via x- To put this strategy into concrete 
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formulation, we proceed first to the momentum space of the partons in the next section, and 
then to the momentum space of the hadrons in the following section. 

3 Correlation between Shower Partons 

Suppose that in a HIC a hard scattering takes place that sends a parton to a momentum 

— * 

k with a large transverse momentum. Such a hard parton generates a shower of partons 
whose momentum- fraction distributions have been determined in Ref. [36]. It should be 
stressed that the shower partons under discussion here are not to be identified with the 
radiated gluons in pQCD. Those shower partons are defined by their recombination with 
each other to form hadrons and their distributions are determined by fitting the appropriate 
non-perturbative FFs in the recombination formalism [36]. The gluons emitted by a hard 
parton, on the other hand, are the products of an evolution process well formulated in pQCD, 
but their hadronization is not well treated because the process is not perturbative. Of course, 
the physics of the shower partons that recombine and the gluons that are perturbatively 
radiated are related in some non-perturbative way, which need not distract us here. Since 
our goal is to determine the correlation between hadrons at intermediate pr, we consider the 
shower partons that are defined for immediate hadronization by recombination in that pt 
range. 

As stated above, those shower partons recombine among themselves to form hadrons 
whose momentum-fraction distributions are the FFs of the initiating hard parton. But when 
the hard parton is produced near the surface of a thermal medium, the shower partons 
can also recombine with the soft thermal partons in the environment to form hadrons at 
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intermediate pr- In the recombination model that is the way how hadron production is 
affected by the medium. Thus for pion production one considers ST recombination, and 
for proton production one has STT and SST recombination, where S denotes the shower- 
parton distribution (SPD) and T the thermal parton distribution. Indeed, it has been shown 
in such an approach that the p/n ratio in Au+Au collisions can be as large as 1 at pr ~ 3 
GeV/c [10]. 

Now to study the correlation between two pions, both in the intermediate pt range, we 
need to consider two shower and two thermal partons in the combination (ST) (ST), where 
the parentheses in (ST) denote the partners that are to recombine. Since only collinear 
partons can recombine, the angle x between the two pions that are formed is the same as 
the angle between the two shower partons, whose momenta magnitudes are generally larger 
than the two softer thermal partons. Let us then concentrate on the two shower partons and 
consider the joint distribution 

{ssy j ' (q u q 2 ) = £E / dkk M k ) ^ (p^>ft) > (f=^'^' a ) } > ( ) 

which is in the form that has been used previously in Ref. [29]. The curly brackets denote 
symmetrization of the momentum fractions of the two shower partons [10], but it is a process 
that is not important in the following, since our emphasis will be on the angular variables. 
For the same reason, we shall not be concerned with £, which denotes the fraction of hard 
partons that emerge from the dense medium to hadronize, and the integral over k, weighted 
by fi(k) that is the distribution of the hard parton % in HIC. ipi and ^2 are the polar angles 
of the j and f shower parton momenta, <fi and q 2 , with reference to the jet axis, k; Pi and 
(3 2 are their corresponding azimuthal angles. Sf ,j are the SPDs. 
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The angle between q\ and q\, denoted by x also, can be expressed similarly as in Eq. (7) 
cos x — cos ipi cos ip2 + sin ipi sin ip 2 cos(/3 2 — Pi) ■ (10) 
If we define, as before, 

^± = ^2±^l, P-=fo-Pl, (11) 

then we have the alternate form similar to Eq. (8), which in turn can be written as 

cosx = A + Bcosp, (12) 

where 

A = i(cos?/>- + c osip + ), (13) 

B — -(costp- — cosip + ), (14) 

and P — /3_ for brevity. Note that A and 5 are linear in costp±, while the corresponding 
terms in Eq. (10) are quadratic in cos^i^- Since x serves as the bridge to the hadron 
momentum space, we want to determine here the distribution in x that correspond to the 
dynamical properties of the jet cone defined with respect to k. 

The joint shower-parton distribution, as expressed in Eq. (9), is based on the assumption 
that the shower partons are dynamically independent, but kinematically constrained by 
momentum conservation 

qi + <?2 < k. (15) 

The constraint on the momentum magnitudes in Eq. (15) does not imply angular constraint 
on ipi and ip2- We shall assume angular independence of the two partons, which implies that 
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their distribution G ? 2('^i, ^2) is factorizable 



G 2 (ii>iM = G 1 (rl> 1 )G 1 (rl> 2 ), 



(16) 



where G\{jp) is the single-parton angular distribution in a jet cone. We assume that it has 
a Gaussian form 



with a width a that is basically unknown, since no reliable theory is calculable at intermedi- 
ate k. We make the assumption here that a is independent of the momentum fraction of the 
shower parton in a jet. That assumption renders significant simplification of our considera- 
tion below, and is a reasonable first step in this exploratory study of autocorrelation. If we 
can relate a to some observable through autocorrelation, we will have achieved the objective 
of probing the microscopic dynamics by phenomenology, at least in the first approximation. 

As stated in the beginning of this section, the shower partons considered in Eq. (16) are 
not directly identifiable as the gluons and their associated low-virtuality partons emitted by 
the hard parton in pQCD because the properties of the partons j and f in Eq. (9) are not 
determined by any perturbative calculation, but by fitting the non-perturbative FFs through 
recombination. Although it is known that a branching process can lead to angular ordering of 
the emitted partons [37, 38], the physics of which has recently been applied to the study of jets 
traversing dense medium [39], the formalism is more reliable at higher jet energy than what 
we are concerned with here at intermediated pr- Besides, hadronization by LPHD is what 
we can avoid when recombination has the advantage of proven phenomenological reliability. 
Note also that we study hadron correlation in jets produced near the surface of a dense 
medium, so the initiating hard partons do not have long path lengths through the medium. 




(17) 
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The only medium effect that we consider is the recombination with thermal partons, which 
does not, in first approximation, influence the angular consideration to follow. For these 
various reasons the angles ipi and ip 2 in Eq. (16) cannot be assigned the angular ordering 
property of the emitted partons in pQCD. For a simple first try we assume factorizability 
in Eq. (16) and the independence of a on the momentum fraction of the shower parton. In 
an extension of this study in the future it would be reasonable to investigate the scenario 
where a increases with decreasing momentum fraction. But for constant a considered here 
the momentum-dependent part factors out and the angular properties are exposed with 
transparency. 

Given the two-parton angular distribution G 2 (ipi, ^2) shown in Eqs. (16) and (17) we can 
calculate the x distribution by allowing for all possible orientations of the two vectors q{ and 
q 2 . Denoting it by H(x), we have 

H(x) = J d cos ip 2 d cos fad/3 G 2 (ip!,ip 2 ) 6 [cos x~ (A + Bcosfi)], (18) 

where Eq. (12) has been used to constrain the four angles ip±,@ and x- Since H(x) is 
invariant under the interchange of q[ and q 2 , we consider only the ranges 

< ipi < ip 2 < a. (19) 

In actual calculation we set a = ir/4. In view of the definitions of A and B given in Eqs. (13) 
and (14), we change the integration variables to ip± and (3. After some algebra we obtain 

H{ X ) = \ r r #- ( c °^-- cos ^^-) _ (20) 

4 Jx Jo [(cos ip- — cos x) (cos x — cos 

where 

i\) m = min(x,2a - (21) 
11 



and 

<7ty + ,lM = exp (-^^)- (22) 

In Fig. 1 (a) we show H(x) for four values of a. To see more clearly the dependence of the 
width of H(x) on the width a of Gi(tp), we normalize the function H(x) by its value at 
X = 0, by defining 

*M = f|. (23) 

It can be shown that in the approximation (j«a, -ff(x) has the limit, as x ~^ 0, -H"(0) ~ 
(7r/4)(7 2 exp(— cr 2 /4). The corresponding normalized distribution H(x) is shown in Fig. 1 
(b). In the a range illustrated, the width of H(x) is ver Y closely given by ou = V^cr. 

It is important to recognize that the factorizable form of (^2(^1,^2) in Eq. (16) does not 
imply an absence of correlation between ipi and ■0 2 . That is because the two shower partons 
with momenta q\ and q 2 are from the same jet, and the angles ip\ and ip2 refer to the same jet 
axis along k. When G 2 (ipi, ^2) is included in the expression for {SS} jj (q u q 2 ) in Eq. (9), the 
contribution to the two-parton distribution p2(<fi, 9*2) is not factorizable. On the other hand, 
p 2 ((fi, ^2) contains a factorizable part Pi(9i)pi(^), in which the two shower partons belong 
to two independent jets that are randomly related to each other. Such a component would 
be cancelled in the calculation of the correlation function ^(gl,^) that can be defined as 
in Eq. (1) for partons. 

Finally, before leaving this section on parton correlation it should be noted that despite 
their appearances the peaks in Fig. 1 should not be confused with the common notion of the 
properties of jet cone. H(x) is the distribution of Xi the angle between q\ and <f 2 , neither 
of which need to be close to the jet axis. That is, x should not be interpreted as the angle 
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between the hard parton and a radiated parton. A collection of all events with a fixed value 
of x is a large set that contains a small subset in which a parton is radiated at an angle x 
with respect to the jet axis. This point is so important that a simple example might be useful 
to make clear the difference. Suppose <fi and q 2 are collinear so that x — 0- However, both 
<fi and q 2 can be directed at some angle ipi = ip 2 that need not be zero, so that x can m 110 
way represent the angle between a radiated parton and the jet axis. Thus H(x) should not 
be taken as a representation of the angular distribution of a radiated parton relative to the 
initiating hard parton. Since x will play a principal role in autocorrelation, this distinction 
underlies the difference between autocorrelation and the correlation based on triggers. 

4 Autocorrelation between Hadrons 

Having obtained the distribution in x, the angle between two shower partons in a jet, we 
can now determine the correlation between two hadrons in a jet. We shall consider only the 
pions, since they are the dominant hadrons in a jet. As we have stated earlier, we consider 
the intermediate pt region where the thermal-shower recombination is most important, i.e., 
3 < p T < 8 GeV/c in central Au+Au collisions [10, 26, 29]. Furthermore, we restrict our 
attention here to only the pions produced at midrapidity so that all momentum vectors 
are nearly transverse to the beam direction. In that case we can simplify our notation by 
denoting the pion transverse momentum px and the parton transverse momentum qx by p 
and q, respectively. The {TS){TS) contribution to the two-particle distribution is then 

dN TSTS 1 r 

— 1 , ,™ , — 3—7— = 7 r-r / dq^q 2 d cos ip^ cos ip 2 d(3 1 d(3 2 T(p 1 - q 1 )T(p 2 - q 2 ) 

p 1 dp 1 dr)id(p 1 p2dp2dr)2d(p2 {P1P2Y J 

{SS}(q u ipi,(3i, q 2 , #2) A(V>i, Pi, 771, 0i; ^2, fo, m, 02), (24) 
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where T is the thermal distribution of exponential form [10], and {SS} is given in Eq. (9). 
The recombination functions, written out explicitly in Ref. [29], have already been integrated 
over, resulting in (a) the momenta of the thermal partons being Pi — qi, and (b) a constraint 
between the pions' angular variables and the shower partons' angular variables ipi,/3i, 
contained in the A function in Eq. (24). We have exhibited only the essential content of 
the (TS)(TS) recombination, where the momenta of the thermal and shower partons that 
recombine are collinear, so all the dynamical characteristics of the problem are in {55}, 
which we have studied in the preceding section, while all the kinematical relationships to 
the observed pions are in A. Although the pseudorapidity variables 771 and 772 appear in Eq. 
(24), it is only a technical step to be taken later to relate their difference to the angular 
differences, 6L and 0_, that appear in the autocorrelation, Eq. (6). Our task now is to first 
relate A{6^,(j)_) to the two-particle distribution in Eq. (24). 

In our study of the autocorrelation our emphasis has been on the dependence on the 
angular differences, Q- and </>_, relegating the px values to the category of other variables 
that are to be integrated over. If we define the full correlation function to be 

C 2 (l,2)=p 2 (l,2)- Pl (l) Pl (2), (25) 

where the arguments symbolize the whole sets of variables of the two detected particles, 
then the LHS of Eq. (24) represents the TSTS component of ^2(1, 2), which we now denote 
by P2 STS (1, 2). Both in the analysis of the experimental data and in our construction of 
autocorrelation, the pt values are integrated over some chosen ranges. In our work here the 
autocorrelation does not depend sensitively on that pr range so long as it is the intermediate 
Pt range where TS recombination is dominant for the formation of a pion. Let us then 
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define the integrated correlation 

r(77i,0i, 772,0 2 ) = J dp 1 dp 2 pip 2 C 2 (l,2). (26) 

With r)± defined by 

V± = V2± Vi (27) 

and (f>± given in Eq. (5), T(i]±, <f>±) receives its contribution mainly from the TSTS compo- 
nent, P2 STS (1,2), expressed in Eq. (24), when i]_ and <p_ become small. 

We now focus on p^ 5T5 (l,2) and, in particular, on the A function in Eq. (24), which 
contains only the angular variables of the parton and pion momenta. As we have emphasized 
earlier, with the autocorrelation A{6_, (f>_) as the aim of our analysis, the angle x is the bridge 
between the parton and pion angular variables that are independent of the coordinate system. 
With that simplification in mind we write A in the form 

A(V>i, (3i, 771, 0i; ip 2 , fa, V2, 02) = 

J dcosx5[cosx - (A + £cos/3_)] 5[cosx - (C + Dcos0_)], (28) 

where A and B are functions of ip±, defined in Eqs. (13) and (14), while C and D are similar 
functions of 6± that follow from Eq. (8), and are explicitly 

C = i(cos#_ + cos#+), (29) 

D = -(cos6_ - cos0 + ). (30) 
2 

It is clear that the two 5 functions in Eq. (28) express the two ways of referring x to the 
two different sets of angles defined with respect to the jet axis on the one hand, and to the 
beam axis on the other. 
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Substituting Eqs. (9) and (28) into Eq. (24), and performing the angular integration as 
in Eq. (18), we obtain H(x) multiplied by some factors involving the momenta magnitudes. 
Such factors are irrelevant to our study, so we normalize it out by defining 

C 2 {v±,<l>±) = v , n t n x ■ 31 
Pending the connection between rj± and 9±, we have from the above process of integrations 
C 2 (9±A±) = J dcos X H(x)S[cosx - (C + Dcoscf)-)}. (32) 

This integration over x need not be performed explicitly, since it is only necessary to de- 
termine the domain of 6± and 0_ that corresponds to each fixed value of x an d denote the 
result as H(9±,<j>J). The final step is to integrate over 9 + as in Eq. (6). Conceptually, it is 
simpler to visualize the average angle 9, defined by 

9 = 9 + /2, (33) 

which should be around n/2, since the two detected pions are restricted to a narrow interval 
at mid-rapidity, corresponding to a jet at roughly n/2 relative to the beam axis. If we denote 
the interval of 9 by 2e, then we have 

A(9 6 ) 1 W 2+e - ~ 

4«-.*-)=4^=&/^'" < *< , *- s --*->- (34 > 

while + is an irrelevant angle in an azimuthally symmetric problem. Experimentally, 
A(9-,(p-) is to be determined by use of Eqs. (6) and (25). 

The transference from 9_ to involves 9 + and is just an algebraic problem. Since 

0i = ^(0 + -0_), 9 2 = l -{9 + + 9_), (36) 
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we can solve for in terms of and 9 + . For 9 not far from tc/2, 9_ behaves essentially 
like rj-. Substituting the exact dependence of 9- on and 9 + into C2(#+,#-,0_) in Eq. 
(34) and then averge over 9, we obtain 



= ^ f'^MH^A-.O)- (37) 

26 Jn/2-t 



This is our main result, except for the 3D display of A(r]_, 0_) after numerical computation. 
The essence of autocorrelation is basically already contained in Fig. 1 (b), where H(x) is 
shown for various values of the jet cone width a. For every value of x there is a set of values 
of rj-, (ft- and 9, whose contour plot gives a representation of A(r)-, 0_). 

In Fig. 2 we plot A(r]_, </>_) for a = 0.2. There is very little sensitivity on the dependence 
on e in Eq. (37); we set it at e = 0.3. The dependence on a is, however, significant, as 
we have already seen in Fig. 1. The general shape of A(r^_,0_) is similar, the peak being 
broader for higher values of a. The utility of this result is, of course, the reverse. When the 
data on A(?7_,0_) become available, we can use A(?7_,0_) to infer what the corresponding 
cone width a should be. To facilitate that deduction, we show in Fig. 3 (a) A(r}-,0) and 
(b) A(0, 0_) for various values of a. We see that the width in is numerically larger 
than that in <f>_. While in the small r/_ region one has ~ 0_, the two variables should 
not be directly compared since r\_ is not an angle. The more important implication of this 
result is that we have in Fig. 3 peaks in the measurable variables and 0_ for four values 
of the theoretical variable a whose magnitude is not known from first principles. Thus any 
experimental information on the autocorrelation can give us direct information on the nature 
of the jet cone. As already mentioned in Sec. I, no trigger-related background subtraction is 
needed for autocorrelation, which is defined in terms of C 2 (l,2). Indeed, experimental data 
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on autocorrelation are already available at low p T [6, 7]. 

Finally, we note that the width of the peak in A(rj-, 0) is wider than that in A(0, 0_), as 
shown in Fig. 3, not because of any flow effect in the longitudinal direction, such as in the 
finding in [18]. We have not considered any medium flow and its interaction with the shower 
partons. Our result follows strictly from a Gaussian distribution for both ipi and ip2 of the 
shower partons relative to the jet axis. The variables 6L and 0_ are differences of angles of 
Pi and f>2 relative to the beam axis, and should not be mistaken to be the angles defined for 
an associated particle with respect to a trigger particle. The example given at the end of 
the preceding section on the angle between parton momenta can be translated now to the 
angle between detected hadrons, thereby illustrating the difference between autocorrelation 
and the correlation between a trigger and its associated particle. 

5 Conclusion 

We have determined the autocorrelation distribution in r\_ and 0_ for two pions produced 
in HIC in the intermediate pt region. We have used the parton recombination model to 
relate the hadronic angles to the partonic angles. By emphasizing thermal-shower parton 
recombination, we have exploited the equivalence of the angle between the two shower par- 
tons and that between the two observed pions. Thus when the data on autocorrelation at 
intermediate px become available, the width of the observed peak can then be related to the 
width of the jet cone, which is a property of the jet physics in HIC that is the goal of this 
study. 

It should be mentioned that our investigation has been focused on the angular relation- 
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ship among the partonic and hadronic momenta. That relationship is independent of pr so 
long as all the transverse momenta in the problem are in the intermediate pr region. Since 
the thermal partons are soft, the shower parton and the pion that is formed by TS recom- 
bination are roughly of the same magnitudes and are collinear. Thus the angle x between 
two shower partons in a jet is the same as the angle x between the two pions, independent 
of their momentum magnitudes. However, we have made use of the simplifying assumption 
that the jet cone has a width a that is independent of the momentum fraction of the shower 
partons. While it is sensible to make that assumption in this first attempt to calculate the 
autocorrelation, we expect the realistic situation to be more complicated. Once that depen- 
dence on momentum fraction is considered, the magnitude of the hard parton momentum 
becomes relevant, and the final result on autocorrelation will exhibit dependence on the pr 
range, even within the intermediate pt region 3 < px < 8 GeV/c. Needless to add, when p T 
goes outside that region TS recombination no longer dominates and the basis for our study 
in this paper will have to be revised. 

Despite the simplifying assumption made in this paper, it will be of great interest to 
compare our result with the forthcoming data on autocorrelation at intermediate px- As far 
as we know, there exists no other theoretical study in the subject that relates the observables 
to the partonic structure in a jet, since no other hadronization scheme has been shown to be 
reliable in the pr range considered. In addition to pions one can also study proton and other 
heavier particles in the jets and their autocorrelations. Even in the pion sector alone it is 
of interest to consider the complications arising from different charge states, since multiple 
shower partons in a jet can exhibit dependence on their flavors. Thus there remains much 
to be learned in the subject, in which the present study is only a beginning. 
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Figure Captions 

Fig. 1. (color online) (a) Distributions in the angle x f° r f° ur values of the width parameter 
a of the jet cone; (b) normalized distributions H(x)- 

Fig. 2. (color online) Normalized autocorrelation function A(?7_,0_) in 3D plot in ?7_-0_ 
for a = 0.2. 

Fig. 3. (color online) (a) Normalized autocorrelation A(r]_,0) for four values of a; (b) 
normalized A(0, 0_). 
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